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Abstract 



Thermoelectric effects through a serial double quantum dot system 
weakly coupled to ferromagnetic leads are analyzed. Formal expressions 
of electrical conductance, thermal conductance, and thermal coefficient 
are obtained by means of Hubbard operators. Results show that although 
thermopower is independent of polarization of leads, figure of merit is re- 
duced by increase of polarization. Influences of temperature and interdot 
tunneling on the figure of merit are also investigated, and observed that 
increase of the interdot tunneling strength results in reduction of the fig- 
ure of merit. Effect of temperature on the thermal conductance is also 
analyzed. 



1 Introduction 

Thermopower through devices fabricated using nanotechnology has been in- 
creasingly studied in recent years [IJ|lliailG3l3IZllEll3linilII]- ReC6nt 
advances in the fabrication of nanodevices and the information showing nan- 
odevices can work as thermal generators with high efficiency |T2l [T3l HI] , have 
increased the importance of such studies. Indeed, violation of the Wiedemann- 
Franz low [T51UB] , ratio of the electrical conductance to the thermal conductance 
is proportional to the operating temperature, in nanodevices results in increase 
of the thermal efficiency and the figure of merit, ZT = s G K vT , where S is 
the thermopower, Gy and T denote the electrical conductance and the tem- 
perature, respectively, k = k c + K p h is composed of the electrical and phonon 
thermal conductance, respectively. Discreteness of energy levels [17], Coulomb 
interactions [T51 Q~5], interference effects [IS], and so on lead to the fact that 
nanodevices have more thermal efficiency than bulky samples. 

In recent years, investigation of the thermoelectrical devices fabricated from 
quantum dots (QDs) has attracted a lot of attention theoretically and experi- 
mentally [fl[2[T![ia[2T][22[23[2![T^ However, thermoelectric effects 
through double quantum dots (DQDs) results in novel phenomena needing more 
studies. Few articles have analyzed thermopower through a DQD in just recent 
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years [27J I2FJ US]- F. Chi and co-workers [25] showed ZT has two huge peaks 
in the vicinity of the electron-hole symmetry points. Trocha and J. Barnas [29] 
reported ZT in DQD systems is enhanced due to Coulomb interactions and 
interference effects. We have recently studied [SU] the effects of Coulomb inter- 
actions and interdot tunneling on the thermopower of a DQD system weakly 
coupled to metallic electrodes. 

In this article, the influences of polarization of the leads, temperature, and 
interdot tunneling on the figure of merit and the thermoelectric conductances of 
a serial DQD weakly coupled to ferromagnetic leads are analyzed. Many-body 
representation introduced in Ref. |31] is used to obtain formal expressions for 
electrical and thermal conductances, thermopower, and figure of merit. Hamil- 
tonian and model used for calculations are described in next section. Section 
3 is devoted to numerical results and in the end, some sentences are given as 
conclusion. 

2 Model and Method 

We consider two single-level quantum dots coupled to ferromagnetic leads. 
Hamiltonian describing the system is given as 

H = ^aka<^ aka C aka + ^ £ i n i<r + X! U i n if n H + U ^ X! n L<y n R<y' (!) 
aha icr i aa' 

+*I>L^ + H -°\ + Y} Vak ° c «k* d ^ + h.c] 

u aka 

where c a fe cr (cJ tfco .) destroys (creates) an electron with spin a, wave vector k, in 
lead a = L,R. Energy levels of the leads are spin-dependent because the leads 
are ferromagnetic. The second and third terms in the above equation describe 
each dot. di cr (d'- )is annihilation (creation) operator of the z'th (i = L,R) dot 
and m a = d ia .dia- Energy levels of each dot are assumed to be degenerate. 
Ui and U12 denote, respectively, on-site and interdot Coulomb repulsions, t 
is the interdot tunneling strength, whereas V a kc denotes the dot-lead coupling 
strength. In order to study the system, the Hubbard operators have been used. 
In this model, the Hamiltonian of the isolated DQD system is diagonalized 
and the creation and annihilation operators are expanded in terms of obtained 
eigenstates. It is obvious that the system has 16 different states shown as 
\N,n > where N — 0...4 is the number of electrons, whereas n denotes n'th 
state of TV-electron configuration. By expanding annihilation operator as 

daa = 2_j (daa^N+l-X-NN+l (2) 
Nnn' 

where the Hubbard operator, X r I !^ N+1 — \Nn >< N + ln'\, describes a tran- 
sition in which an electron inside the DQD system is annihilated, the whole 
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Hamiltonian is given as 

H = ^ £akgC i aka C aka + S ^ j E Nn h T ^+ ^ [KiA : cr((i a(T )^r I 5 v+1 C^ fe(T X^ r '] v+1 +i/.C] 
afccr JVji akaNnn' 

(3) 

where /i^ = |JVn >< Nn\. 

Now, population number Pjvn, the probability of being in the state \Nn >, 
is computed by means of density matrix approach. Coupling to the leads is 
considered so weak that non-diagonal elements of the density matrix are ignored 
because they are proportional to V£ kcr . Using Markov approximation and wide 
band limit, the time evaluation of the population numbers is given as |31] 

°- = X^[ _r roi>^|i»>- P ° 1 + r fin>^|oi> i:, i«] ( 4a ) 

an 

( jl ^~]~\F?Nn>->\N-ln'> ^UVn>^-liV+ln'>]-^V""^ (^-0 
an' 

r|jv_i„'>_ > .|iv n >-PAr-ln / + F" N+lnl> ^ Nn> P N+ln > 

X^ r f41>^|3n>- P 41 + r f3n>^|41>^3n] (4c) 



dt 

an 

dP Nn 



dt 

an 

It is obvious that there is one configuration for zero- and four-electron state. 
The transition rates are 

^"\Nn>^\N+ln'> = I (^aa )nn+1 P/a i^N+ln' ~ E Nn ) (5a) 

a 

^\k+\n'>-t\Nn> — t 5Zr^|(d Qcr )^ v+1 | 2 / Q (E N+ln > — E Nn ) (5b) 

a 

where T° — 2ir J^kea Waka | 2 is the spin-dependent tunneling rate, f a (x) — 
(1 + exp((x — M Q )/fcr Q )) _1 is the Fermi-Dirac distribution function in which 
\i a and T a stand for the chemical potential and the operating temperature of 
the lead a, respectively, and f~ = l — f a . The charge and energy currents are 
computed by solving the Eqs. (|4]) in the steady state situation ( J[ n = 0) 

I a = -e ^Nn^lN^P^sgniN' - N) (6a) 

N,N',n,n' 

Q a = E T \Nn>^\NW>( E \N'n>> - E \Nn>)P\Nn>Sgn(N' - N) (6b) 
N,N',n,n' 

where sgn(x) is a signum function. 

To compute the thermoelectrical characteristics of the system, the linear 
response regime is used. We assume that the left lead is slightly hotter than the 
right one (T& —Tr + AT), and = — eAV, so that the charge and energy 
currents are given as follows 

I a = G V AV + G T AT (7a) 
Q a = MAV + KAT (7b) 
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where Gy is the electrical conductance and Gt is the thermal coefficient. Ther- 
mopower is defined as minus ratio of induced voltage to applied temperature 
gradient under condition that the charge current is zero, so we have S = 
Putting I L = l/2(I L — I R ) and expanding the fermi-Dirac distribution function 
as f L ( x ) = f R (x)-x/Tf'(x)AT + eAVf(x) where f'(x) = df(x)/dx, one can 
easily obtain: 

e 2 

Gv = — 2^ z2i a2 p N7i^^[\(d La ) , pP N+1 \ 2 f'(E N+ln , - E Nn )+ (8a) 

Nnn' o 

|(dZff)jV— Uvl /'(-^iVn — £-JV-ln')] 

Gt = y^Pj V „r^[|(d L(T )^ +1 | 2 ( ' N+ln — f'(E N+ i n > - E Nn ) + 

Nnn' cr 



<„ . 2 (E Nn - E N -i n i) , 



(8b) 



|(d£ff)j\r-livl y / ( E Nn - ^JV-ln')] 



Using above equations and Onsager relation, the thermal conductance is com- 
puted as [T7] 

k c = [if - S 2 G V T] (9) 

where 

^ = -^z X! X] jJ 'V" r ^[l(^^)jVjv+il 2 ^ jV+1 " r ^—f{E N+w - E Nn )+ 

Nnn' o 



if j \n'n \2 ^Nn E N _ ln >) 2 , 



(10) 



For simulation purpose, we assume that = T R = To and spin-down 
tunneling rate is equal to = T R = oFo where < a < 1. It is obvious that 
a = 1 denotes normal metallic electrodes whereas a = stands for half-metal 
leads, we also set n p h — 3kq where kq = 3h B T is the quantum of thermal 
conductance [35] and, assume that the single electron levels in the QDs are 
degenerate. In recent years, the phonon contribution in the transport through 
QDs has been extensively studied experimentally and theoretically |33j [3H [35] . 



3 Results and discussions 

Figure of merit (ZT) as a function of the QDs' energy level and the temperature 
is plotted in fig. 1 for different as. Similar plot was previously presented in 
Ref . [18] about a QD coupled to Ferromagnetic electrodes or in Ref. [29] about a 
DQD coupled to external electrodes. One can observe that increasing a results 
in the reduction of ZT. On the other hand, ZT has some peaks whose intensities 
decrease with increase of temperature, but become wider. The results also show 
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that ZT approaches zero in — 3 < Si < and high temperature (T > 4k). With 
respect to the fact that ZT is a function of the Seebeck coefficient, the electrical 
conductance, and the thermal conductance, the evaluation of them will help us 
realize the behavior of ZT. 

Figure. 2a describes the electrical conductance as a function of the energy 
level in different situations. From eq. (8a) it is obvious that Gy has a peak 
whenever the necessary energy for transition from TV-electron state to N ± 
1-electron state equals to the Fermi energy of the leads (fi a =0). In low 
temperature (T — 2k), f'(x) is a Lorentzian function with narrow width, thus 
the peaks and the Coulomb valleys are clearly observed. There are four peaks 
because four electron can exist in the system. With increase of the temperature, 
the intensity of the peaks reduces whereas they become wider. Furthermore, the 
Coulomb valleys are vanished because increase of the temperature weakens the 
Coulomb blockade effect. It is straightforward to show that Gy = -MpGo where 
G is the electrical conductance of a double quantum dot coupled to normal 
metal electrodes. Therefore, the polarized electrodes decrease the electrical 
conductance of the system. The thermal coefficient as a function of energy is 
shown in fig. 2b. One can easily observe from eq.(8b), for energies in which Gy 
has a maximum Gt and as result, S become zero because of Enu ~ £W±in'- 
In these energy points, although the transition from iV-electron state to N ± 
I-electron state results in the electrical current, no net energy transports in 
the process. Indeed, the temperature gradient does not have any significant 
role in producing the charge current in resonance energies. In addition, the 
thermal coefficient becomes zero in some other energies. From eq.(8b), it is 
obvious when the system is in a certain state {Pnu — 1) and E^+in' — E^ n = 
— (Enu — EN-in>), then Gt = 0. These energy points are so-called electron- 
hole symmetry points and the effect has been recently studied in a multilevel 
QD [36] . In symmetry points, the current transported from the lead to the QD is 
equal to the current transported in opposite direction, so that the thermopower 
becomes zero. Indeed, both electrons and holes participate in producing the 
current with different signs. It is worth noting that Gt goes from positive to 
negative values in resonance points whereas it behaves reversely in symmetry 
points. This effect, bipolar effect, has been recently reported in a multilevel 
quantum dot [36]. It comes from the fact that in each side of it, a kind of charge 
carrier (electron or hole) participates in transport. This leads to the oscillation 
of the thermopower. Increasing temperature gives rise to widening Gt but its 
intensity will be reduced. The influence of the polarization of the leads on 
the thermal coefficient is the same as Gy and as a result, the thermopower 
is independent of the a. This result was previously reported about a single 
level QD coupled to ferromagnetic leads [20 . However, results presented in 
Ref. |18j show that the thermopower weakly depends on the polarization of the 
leads. This difference comes from the different models used for studying the 
system. Indeed, in nonequilibrium Green's function formalism used in Ref. |18j . 
the broadening of the QD level due to coupling between the QD and the leads 
is taken into account. That leads to the dependence of the thermopower on 
the magnetic polarization. In rate equations approach, the broadening of the 
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QD level due to coupling is disregarded. Figure of merit is zero in resonance 
and symmetry pints as one can see in fig. 1 due to ZT oc S. Furthermore, ZT 
decreases with increase of temperature because increasing temperature results 
in reducing S and Gy. On the other hand, although the thermopower is spin- 
independent, the figure of merit is reduced by increase of polarization because 
of ZTocGy. 

The thermal conductance as a function of temperature and energy level is 
shown in fig. 3. Results show that the thermal conductance takes the highest 
values in electron-hole symmetry points when the temperature is sufficiently 
high. In these points, electrons and holes participate in the transfer of charge 
and heat with the same weight. Because the current carried by holes is in the 
opposite direction with one carried by electrons, the charge current becomes 
zero in these energies so that one can observe as valleys of Gy in fig. 2a. 
In contrast with the charge current, electrons and holes carry the energy in 
the same direction resulting in the appearance of the peaks in k c . As expect, 
increasing temperature gives rise to the increase of k c . In low temperature, 
the peaks of Gy and happened in nearly same energies. Indeed, the 

significant difference between Gy and n c is observed in high temperature. In 
low temperature, when the transition energy is equal to resonant energy, one 
peak develops in both Gy and K c , resulting due to the transfer of current and 
heat through the system by electrons jTS] . 

Figure. 4 describes figure of merit as a function of the polarization of the 
leads and the interdot tunneling strength. It is observed that increase of the 
polarization of the leads results in the reduction of ZT. Results also show that 
the interdot tunneling leads to decrease of ZT. Indeed, From fig. 4 one can 
approximately consider ZT oc at so that ZT obtained from strong a (weak 
polarization) and weak t is equal to result obtained from weak a (strong po- 
larization) and strong t. It is important to note that the eigenvalues of the 
DQD are a function of the interdot tunneling, so the position of symmetry and 
resonance points depends on the interdot tunneling strength. It is worth noting 
that increase of interdot tunneling results in reduction of Gy [30] . Furthermore, 
although the thermopower is spin- independent, the electrical conductance is re- 
duced by increasing polarization, and as a result the figure of merit reduces when 
the polarization of leads increases. In addition, intra- and interdot Coulomb re- 
pulsions affect the figure of merit significantly. Recently, we have shown that 
the increase of Coulomb repulsion results in the increase of figure of merit [30] . 
Indeed, the level spacing is increased by increase of Coulomb repulsion and as 
a result, the bipolar effect will be reduced. The same effect was previously 
reported about a multilevel QD [36] . 

4 Conclusion 

In this paper, the thermopower in a serial double quantum dot system weakly 
coupled to ferromagnetic leads is studied using Hubbard operators. Formal 
expressions for thermal conductance, electrical conductance, and thermal coef- 
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ficient are obtained using density matrix approach. The effect of temperature, 
interdot tunneling, and polarization of leads on the thermoelectrical character- 
istics of the system is examined. Results show that increase of temperature 
results in decrease of electrical conductance intensity whereas it becomes wider. 
It is found that although the thermopower is independent of polarization, the 
figure of merit is reduced by increase of polarization. Influence of the interdot 
tunneling on the figure of merit is also analyzed. 
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Figure captions 

Figure. 1: Color map of figure of merit. Parameters are: t = 0.2meV, 
Ui = 2meV, U 12 = ImeV and T = lOfieV. 

Figure. 2: Gy and Gt as a function of energy. Solid (T = 2K, a = 1), 
dashed (T = 5K,a = 1), dash-dotted (T = 2K,a = 0.5), and dotted (T = 
5K, a = 0.5). Other parameters are the same as fig. 1 . 

Figure. 3: Color map of thermal conductance. Parameters are the same as 
fig. 1. 

Figure. 4: Figure of merit as a function of polarization and interdot tunnel- 
ing. Parameters are the same as fig. 1 except T — 2K and Si = —1.5meV. 



9 






x 10 
1.5 



0.5 



I — u 



0.5 



2 4 6 



T(k) 



2 4 6 

T(K) 



Figure.l 



10 




11 




12 



8 




t(meV) 

FigureA 



13 



